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ABSTRACT 

The gas-phase metallicity distribution has been analyzed for the hot atmospheres of 
29 galaxy clusters using Chandra X-ray Observatory observations. All host brightest 
cluster galaxies (BCGs) with X-ray cavity systems produced by radio AGN. We find 
high elemental abundances projected preferentially along the cavities of 16 clusters. 
The metal-rich plasma was apparently lifted out of the BCGs with the rising X-ray 
cavities (bubbles) to altitudes between twenty and several hundred kiloparsecs. A 
relationship between the maximum projected altitude of the uplifted gas (the “iron 
radius”) and jet power is found with the form i?Fe oc The estimated outflow rates 

are typically tens of solar masses per year but exceed 100 Mq yr“^ in the most powerful 
AGN. The outflow rates are 10% to 20% of the cooling rates, and thus alone are unable 
to offset a cooling inflow. Nevertheless, hot outflows effectively redistribute the cooling 
gas and may play a significant role at regulating star formation and AGN activity in 
BCGs and presumably in giant elliptical galaxies. The metallicity distribution overall 
can be complex, perhaps due to metal rich gas returning in circulation flows or being 
blown around in the hot atmospheres. Roughly 15% of the work done by the cavities 
is expended lifting the metal-enriched gas, implying their nuclear black holes have 
increased in mass by at least ^ 10^ Mq to 10® Mq. Finally, we show that hot outflows 
can account for the broad, gas-phase metallicity distribution compared to the stellar 
light profiles of BCGs, and we consider a possible connection between hot outflows 
and cold molecular gas flows discovered in recent ALMA observations. 
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1 INTRODUCTION 

Radio active galactic nuclei (AGN) drive shock waves and 
inflate giant bubbles into the hot, X-ray atmospheres per¬ 
vading galaxies and clusters (for reviews see McNamara & 
|Nulsen|[2007 2012 Fabian|[2012 1. The enthalpy released as 
the bubbles rise through hot atmospheres is apparently able 


Grand! fc Molendi|2001 1 . The hot gas became enriched ap¬ 


parently during the early development of clusters, primarily 
through ejecta from core-collapse supernovae (SNe II) as¬ 
sociated with star formation. The atmosphere immediately 
surrounding the brightest cluster galaxy (BGG) is chem¬ 
ically enriched above the cluster mean, approaching and 


regulate the rates at which hot atmospheres condense into 
stars I Bower et al.|2006 Sijacki fc SpringeH2006 Croton et 


to offset radiative cooling (Birzan et al. 2004 Vernaleo & 

Dupke & White||2000| lUe Grandi & MolendipOOlMTamura 

Reynolds 2007 Mendygral, O’Neill, & Jones 20111 and to 

et al.|2001 

1 . This plasma is apparently stripped debris from 


al.|2006 1. The rising X-ray bubbles or cavities also drag cool, 
metal-rich, plasma away from the host galaxy and deposit 
it on larger scales (Kirkpatrick, McNamara, & Cavagnolo 
2011). The rate with which this hot, metal-rich gas flows is 
largely unexplored, and thus may be a significant aspect of 
radio-mode feedback. 


orbiting galaxies and stellar ejecta from the BCG itself. Be¬ 
ing the dominant stellar structure in the core of a cluster, a 
BCG enriches cluster cores through stellar winds and type la 
supernovae (SN la), a process generating greater quantities 
of iron compared to SN II (De Grandi et al.|2004 Tamura et 


ar]|2004 de Plaa et al.|[2006 K A dichotomy therefore exists 

between the a-element enhanced outskirts and the iron-rich 


Cluster atmospheres are replete in heavy elements with 
average metallicities of approximately 1 /3 of the solar value 


iMushotzky et al.|1996 Mushotzky fc Loewenstein|1997 De 


Were the iron-rich ejecta produced by a BCG lying 
in a static atmosphere, the stellar light and atmospheric 
metal concentrations would track each other. However, sev- 
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2 C. C. Kirkpatrick & B. R. McNamara 


eral studies have shown that the radial distribution of iron 
and other metals in cool core clusters is more broadly dis¬ 
tributed than the stellar light profile iRebusco et al.|[2005 


2006 David &; Nulsen||2008 Rasera et al.||2008 l. Evidently, 

the metals produced by BCGs are transported to higher al¬ 
titudes. Likely agents include bulk and turbulent motions 
driven by halo mergers and active galactic nuclei (AGN) 


(|Gopal-Krishna & Wiita 

2001 

2003 

Briiggen & Kaiser 

2002 

Omma et al. 2004| [Mol 

et a 

.[[20071 [Roediger et al.[ 

2007 


Pope et al.pOlo'l. A growing body of data implicate AGN. 


The Virgo and Hydra A clusters are archetypes. In M87, 
bright, X-ray filaments that are metal enriched with re¬ 
spect to the surrounding gas are draped along the radio 
jets, indicating that the filaments were lifted out by the jets 


dra A but on vastly larger scales exceeding 100 kpc ([Nulsen 

et al.|2002[ 

Simionescu et al.|2009[ [Kirkpatrick et al.|2009a[ 

Gitti et al. 

2011 

I. Hydra A is experiencing a much more 


powerful AGN outburst compared to the current outburst 
in M87, which is why the metal rich gas is seen at such 
large altitudes. This phenomenon is observed in rich clus- 


ters (Sanders et aL|2005 

Doria et al.|2012 

) and 

groups, e.g.. 

HCG 62 (Morita et al. 

2006 Gu et al. 

2007 

Gitti et al. 

2010 Rafferty et al. [[20131. A study of ten clusters experi- 


encing radio-AGN activity (Kirkpatrick et al. 20111 revealed 


a power law relationship between jet power and the largest 
projected radius to which metal-enriched gas is observed. 
Indications are that persistent AGN activity may be largely 
responsible for enriching the cluster core. 

Here we present an analysis of 29 galaxy clusters drawn 
from the Chandra X-ray Observatory archive. The sample 
was chosen based on the level of AGN activity exhibited by 
X-ray cavities and availability of long exposures needed for 
accurate gas phase metallicity measurements. In section 2 
and 3 we present the entire sample of clusters and discuss 
preparation of the X-ray data. Section 4 is our analysis and 
discussion of the sample. Our main goals are to further re¬ 


fine the relationship established by Kirkpatrick et al. (20111 


between jet power and the extent metals are uplifted out of 
the cluster centre and to test whether or not an AGN can be 
the driving mechanism behind uplifting all the gas originally 
associated with the BGG and creating the broad abundance 
peaks we currently observe. Sections 5 is the summary of our 
results. Throughout this paper we assume a AGDM cosmol¬ 
ogy with Ho = 70 km s“^ Mpc”'^, Dm = 0.3, and Da = 0.7. 
All uncertainties are quoted at the 67% confidence level. 


2 CLUSTER SAMPLE 


Twenty nine clusters, fo und in Table [Tj w ere selected from 
[Rafferty et al.| ( |2006[ ) and [Dunn fc Fab^n| ( |200^ , most with 
clear AGN activity in the form of X-ray cavities. The sam¬ 
ple is divided into two subsamples, dubbed “High Qual¬ 
ity” (HQ) and “Extended” (EX). The High Quality sample 
alone is used to comparing jet power to uplifted metals. This 
sample is composed of clusters with Chandra images deep 
enough to provide accurate gas metallicity measurements 
projected on scales smaller than the radio/cavity system. 
It includes the original 10 clusters with which the trend 
between the jet power and iron radius was found pSF] 
Patrick et al.|[20TT|, and an additional six with deep Chan¬ 


dra data that have since become available. The “Extended 
Sample” includes 13 additional systems most of which have 
X-ray cavity jet power measurements, and central metallic¬ 
ity peaks, but whose Chandra images are too shallow to map 
the metallicity distribution on fine scales. The two samples 
are combined to analyze global metallicity properties and 
for the analysis of growth rates of super-massive black holes 
(SMBHs) in BCGs, discussed in Section 4.4. 


3 DATA PREPARATION 

The X-ray data were taken with the Chandra X-ray Ob¬ 
servatory using either the ACIS-S3 back-illuminated CCD 
or the ACIS-I front-illuminated CCD array. Data prepara¬ 
tion was carried out using version 4.1.2 of both Cl AO and 
the calibration database. When present, background flares 
were eliminated using the lc_CLEAN routine. Over the full 
energy range (0.3-10 keV), a light curve was extracted from 
the level 2 event file. Time intervals with photon count rates 
varying greater than 3-sigma from the mean were removed. 
The corrected exposure time for the combined observations 
can be found in Table Time-dependent gain and charge 
transfer inefficiency corrections were applied. Background 
subtraction was carried out using blank-sky background files 
normalized to the source count rate in the 9.5 — 12 keV en¬ 
ergy band. Coordinates were reprojected to match the new 
background file to the source image in order for background 
spectra to be extracted from the equivalent spatial regions. 
Targets with multiple observations were added together to 
properly identify faint point sources using wavedetect. Each 
point source was confirmed manually and removed from the 
final level 2 event file. 


3.1 X-ray Spectral Modeling 


Most observations were taken at a focal plane temperature of 
— 120°C, allowing us to create weighted response files using 
the CIAO tools mkacisrmf and mkwarf. For observations 
made before 2000 January 29, mkrmfwas used instead. We 
modeled all spectra using a single temperature plasma model 
plus photoelectric absorption unless otherwise noted. Using 
a WABSxMEKAL model over the energy range 0.5-7 keV, 
temperature, abundance, and normalization are allowed to 
vary. The column density was set to the values quoted by 
[Dickey fc Lockmanj ( [1990[ ) . Abundance ratios (the fraction 
of an element of the gas compared to its fraction in the 


sun by number) were set as in Grevesse & Sauval (19981, 
based on solar photospheric data. These ratios are tied to 
the abundance of iron, as it is the prominent feature in the 
X-ray spectrum. 

In cases where single temperature models provided an 
inadequate fit to the spectrum, multiple temperature model 
fits were performed eg., ( Kirkpatrick et al.|2009b l. The fail¬ 
ure of the single temperature model at the centre, where 
the temperature structure becomes complex, usually reveals 
itself as a sharp decrease in central metallicity instead of 
the usual rising metallicity peak (|Buote|[2000[ [Molendi &[ 


Gastaldello 20011. The single temperature model tends to 


underestimate emission at energies above 4 keV. This is¬ 
sue has been noted in the cores of clusters and in some 
AGN outflow regions. An example of a typical spectrum 
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Table 1. The Sample. 


System 

Sub-sample 

Exp. Time 
(ks) 

z 

p a 
-^jet 

(lO"'^ erg s-l) 

Mb'" 

Re'’ 

(kpc) 

rp d 

J- core 

A133 

HQ 

103.6 

0.0566 

620 

-22.64 

25.0 

IT 

A262 

HQ 

139.4 

0.016 

9.7 

-21.07 

12.3 

2T 

Perseus 

HQ 

957 

0.018 

150 

-22.65 

— 

IT 

2A 0335 

HQ 

102.4 

0.035 

24 

-23.04 

— 

2T 

A478 

HQ 

52.4 

0.088 

100 

-23.54 

— 

IT 

MS 0735 

HQ 

476.5 

0.22 

6900 

-22.18 

— 

IT 

Hydra A 

HQ 

181.9 

0.055 

430 

-23.05 

39.3 

IT 

HCG 62 

HQ 

167 

0.014 

3.9 

-19.97 

— 

IT 

A1795 

HQ 

205.8 

0.063 

160 

-22.59 

— 

IT 

A1835 

HQ 

193.7 

0.253 

1800 

-23.27 

— 

IT 

A2029 

HQ 

97.7 

0.077 

87 

-19.87 

25.1 

IT 

A2052 

HQ 

617.3 

0.035 

150 

-21.74 

14.1 

2T 

A2199 

HQ 

156.5 

0.030 

270 

-23.01 

14.2 

IT 

Sersic 159 

HQ 

107.6 

0.058 

780 

-22.34 

— 

IT 

A2597 

HQ 

134.4 

0.085 

67 

-22.51 

— 

2T 

M87 

HQ 

480.1 

0.0043 

6 

-22.05 

— 

2T 

A85 

EX 

38.4 

0.055 

37 

-22.76 

31.3 

IT 

PKS 0745 

EX 

40.2 

0.103 

1700 

-20.53 

— 

IT 

4C 55.16 

EX 

95.8 

0.241 

420 

-22.14 

— 

IT 

RBS 797 

EX 

49.7 

0.35 

1200 

— 

— 

IT 

A2390 

EX 

110.6 

0.23 

* 

— 

— 

IT 

Zw 2701 

EX 

122.5 

0.21 

* 

-22.25 

— 

IT 

Zw 3146 

EX 

83.3 

0.29 

* 

-22.81 

— 

IT 

MACS J1423 

EX 

133.6 

0.54 

1400 

— 

— 

IT 

MKW 3S 

EX 

57.3 

0.045 

410 

-22.09 

18.8 

IT 

Hercules A 

EX 

96.8 

0.15 

310 

-21.51 

— 

IT 

3C 388 

EX 

38.4 

0.091 

200 

-22.98 

— 

IT 

Cygnus A 

EX 

225.8 

0.056 

1300 

-21.83 

— 

IT 

A4059 

EX 

91.9 

0.048 

96 

-22.88 

23.5 

IT 


“ Jet power from [Rafferty et al.| l |2006l l 
^ Absolute magnitudes derived from HyperLeda catalogue 
Effective radius from HyperLeda catalogue 
Temperature model used for X-ray core 


is presented in the left panel of Figure The residuals 
show data above 4 keV systematically above the single 
temperature model. When present, the effect was identified 
and corrected by adding a second temperature component 
(WABSx (MEKAL+MEKAL)) for the metallicity maps and 
profile fits. The second component’s temperature and nor¬ 
malization were allowed to vary with the abundance tied be¬ 
tween the two component. The two-temperature fit shown 
in the right panel of Figure is shown for example. We 
then applied the two-temperature model (2T) to all spectra, 
including those with acceptable single temperature model 
fits. In those cases we found the second temperature compo¬ 
nent consistent with zero, or an improvement in chi-squared 
statistics that, when performing an F-test, was shown to be 
insignificant. We list in Tablewhich temperature model is 
used for each cluster. When needed, the 2T model is only 
used in the inner 1 or 2 bins of the profile or central region 
of the map, except in M87, where the entire profile along 
the jet region uses the 2T model. 


3.2 Systematic Errors 


To control for systematic errors, we explored the effects of 
misestimation of the foreground hydrogen column density on 
our metal abundance measurements. We initially assumed a 
uniform neutral hydrogen column density across the field 
of view. However, this assumption is invalid toward fields 


projected onto Galactic molecular clouds (Bregman et al. 


2003 Bourdin & Mazzotta 20081. The spectrum presented in 
the left panel of Figure [2| is an example of a poor spectral fit 
owing to variations in column density across the field of view. 
Based on images from the IRAS all sky survey, we found 
that clusters with unacceptably high residuals below 1 keV 
were all projected toward high Galactic column densities. 
The right panel of Figure represents the improvement in 
spectral fitting when column density allowed to adjust itself 
across the field of view. 


4 ANALYSIS & DISCUSSION 
4.1 Metallicity Maps 

We present projected metallicity maps for the High Quality 
sample of clusters in Figurej^ The spectra were binned using 


© 0000 RAS, MNRAS 000, 000-000 












4 C. C. Kirkpatrick & B. R. McNamara 




Energy [keV] Energy [keV] 

Figure 1. Left: A typical, multi-observation X-ray spectrum poorly fit by a single-temperature model. Right: The same spectrum is 
better fit at energies above 3 keV when a second temperature component is added. 



0.5 1 2 5 

Energy [keV] 



Energy [keV] 


Figure 2. Left: A typical, multi-observation X-ray spectrum, where the black and red data are from ACIS-S and ACIS-I, respectively. 
It is poorly fit when the average column density is underestimated. Right: Allowing the column density to vary when large molecular 
clouds are present along the line of sight improves the fit at energies below 1 keV. 


a weighted Voronoi tessellation (WVT) algorithm ( [Cappel-] 
lari & Copin 2003 Diehl fc Statler|2006 1. We began by defin¬ 


ing a minimum signal-to-noise (S/N) of 100 per bin. The 
value was increased to between 150 and 200 when possible. 
We optimized the analysis to achieve bin sizes that minimize 
uncertainties but provide sufficient sampling to recover spa¬ 


tial trends in the metallicity distribution on angular scales 
below that of the radio/cavity system. 

The WVT algorithm searches for the highest S/N pixel 
in the input image, which usually lies near the cluster centre. 
The bin is grown by merging with the highest S/N neighbour 
and the centroid is then recalculated. The process repeats 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 



0.4 0.5 0.6 0.7 0.8 0.9 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 


Figure 3. Metallicity maps for the high quality sample. The green sectors overlaid are along the axis of the launched X-ray cavities. 
These bins are where the spectra are extracted when creating the on-jet metallicity profile. 


until the bin reaches the target S/N. Bins are created sim¬ 
ilarly until the cluster is mapped at the desired S/N per 
pixel. Maps of the more shallowly exposed objects from the 
Extended Sample were generated at constant spatial reso¬ 
lution rather than a targeted S/N. We chose bin sizes to 
insure more than a single bin per quadrant on the sky and 
S/N ranging between 45 and 75. 

We do not use the metallicity maps to measure the iron 
radius discussed below. The maps in Figure]^ are intended 
to show how the gas phase metallicity changes with altitude 
and azimuthal angle. In each panel of Figurewe have indi¬ 
cated the spectra extraction regions for abundance profiles 
along the axis of the cavities system. The orientation and 


width of the bins approximates the orientation and width of 
the cavity systems. 


Inspection of the maps reveal plums of high metallic¬ 
ity gas extending from the centre. Striking examples include 
Hydra A, Abell 1795, Abell 478, and particularly MS0735, 
which reveals a prominent “rust streak” to an altitude ap¬ 
proaching 300 kpc along the cavities. Lacking direct veloc¬ 
ity measurements, we assume the metals track outflowing 
plasma lifted by the bubbles. This assumption may be tested 


in future with the combination of Astro-H (Takahashi et al. 


20121 and Athena. 


Note that high metallicity gas is found preferentially, 
but not exclusively, along the bubbles. High metallicity 
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Figure 3. - continued 


patches are found away from the bubbles as well. A good 
example is Hydra A, which shows a high metallicity plume 
away from the cavities associated with an X-ray filament 
( Kirkpatrick et ^|2009a f These filaments may be gas re¬ 
turning in a circulation flow ( [Brighenti fc Mathews|2006[ ) or 
gas being blown around in the dynamic intracluster medium 


( Heinz et al.||2006 Morsony et al. 2010 Ruszkowski & Oh 
|2010[ [Simionescu et al.|2012[ ). 


The situation in the Extended Sample, for which we 
have not shown the metallicity maps, is less clear. The spa¬ 
tial bins are too large to reveal any convincing evidence for 
the variation of metallicity along the cavity axis. Due to 
these large uncertainties in individual systems, they are not 
included in the trend analysis below. 


4.2 Updated Pjet - Rpe Scaling Relation 

With evidence from the metallicity maps for metal-rich 
plumes, we have estimated the maximum altitude of the 
metallicity enhancement in the ICM. The existence and size 
of the enhancement was found by comparing the metallicity 
profiles of each cluster along and orthogonal to the X-ray 
cavities. The profiles were created using semi-annular bins 
with opening angles lying between 45 and 90 degrees for 
most clusters. One profile enclosed the region along the X- 
ray cavities and/or where extended radio emission from a jet 
is found (“on-jet”). The second profile in the orthogonal di¬ 
rection samples the undisturbed atmosphere (“off-jet”). We 
assume the metallicity profile of the undisturbed atmosphere 
represents the average prior to recent AGN activity. 
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To find deviations between the on-jet and off-jet pro¬ 
files, each semi-annular bin included enough counts to re¬ 
duce the measurement uncertainties to approximately 10%. 
This criterion was crucial to preserve spatial resolution be¬ 
cause larger radial bins tend to average away small scale 
metallicity fluctuations due to the overall metallicity gradi¬ 
ent. This criterion was met by extracting profiles with an 
average of 9 bins with minimum S/N of 140 for each bin. 

Profiles for the High Quality sample are found in Fig¬ 
ure]^ The on-jet profiles are represented by the circular data 
points. The off-jet profiles are represented by triangles. The 
dotted vertical line represents the location of the “iron ra¬ 
dius”. We define this point by counting bins outward from 
the centre until we reach the bin where the la error bars 
do not overlap and the regions beyond this radius are ei¬ 
ther indistinguishable between sectors or the off-jet region 
is higher in metallicity. 

With a smaller subset of data, a relationship has been 
found showing the iron radius scales with the AGN jet power 
( Kirkpatrick et al.|2011 1. We refined this trend using the en¬ 
tire high quality sample which we expanded here to include 
additional low-power and high-power systems. The new jet 
power vs. iron radius plot is found in Figure[^ The jet pow¬ 
ers used here were taken from cavity measurements by |Raf-| 
ferty et al. (20061, except for Sersic 159/03, which were re¬ 


calculated by including new data from a longer exposure 
time observation (ObsID 11758.) The “error bars” for Rpe 
represent the width of the bin in the metallicity profiles (Fig¬ 
ure]^. After performing a least-squares linear regression, the 
relationship is found to be. 


i?Fe = (62 ± 26) X P^«±0.06) 


( 1 ) 


jet 

where jet power is in units of 10^"^ erg s“^. The standard 
deviation is approximately 0.87 dex. Increasing the sample 
size from 10 to 16 changed the relationship imperceptibly. 
The predicted iron radius is slightly higher, but within the 
error of the sample the new relationship is indistinguishable 


from that in Kirkpatrick et al. (20111. 


We have made no attempt to correct Figure for pro¬ 
jection. Our pV calculations assume the outbursts lie in the 
plane of the sky so that the line-of-sight diameter of a cav¬ 
ity is similar to the projected major and minor axis. It is 
important to understand what effect this could have on our 
results and final conclusion. 

The iron radii suffer similar uncertainties to the pressure 
and buoyancy timescales. The buoyancy time is defined as 
the time for a cavity to rise at its terminal velocity, 

tbuoy = Ry/SCl2gV, (2) 

where R is the cavities projected distance from the cluster 
centre, S is the cross section of the cavity, C is the drag co¬ 
efficient, and V is the cavity volume ( [Churazov et al.||2000[ 
2001). In general, cavities are oriented at oblique angles, 


which underestimates the altitude and buoyancy time, and 
overestimates pressure. In the extreme, metal asymmetries 
in systems with jets aligned near the line of sight could be 
foreshortened enough to avoid detection. On average, the 
cavities are expected to be oriented 45 degrees or more 
because cavities projected directly along the line of sight 
against the nucleus are difficult to detect. We would then 
expect the cavity altitude and iron radius to be underesti¬ 
mated by 30%. Any attempt to correct for projection would 


shift the Pjet — Rpe scaling relation to larger distances, but 
would otherwise leave our conclusions unchanged. In Fig¬ 
ure]^ most points on the trend would shift to higher powers 
and a greater iron radius. The overall effect is it may take 
less power to move material to further distances than our 
relationship implies. 

As discussed in (Kirkpatrick et al. 2011| , qualitative 
agreement was found between our observed scaling relation 


and jet simulations. Morsony et al. (2010) studied the spa¬ 


tial impact jets have on their environment, finding a de¬ 
pendence of i? oc a shallower relationship than what 


we find in this work. In the work of Gaspari et al. (2011), 
with a simulated jet power of 5.7 x erg s“^, the gas is 
lifted to 200 kpc from the cluster centre. This is larger than 
what our scaling relationship predicts by about la, but the 
overall agreement is encouraging. We expect that attempts 
to match our scaling using simulations will levee restrictive 
constraints on jet physics. 


4.3 Correlations with Cavity Enthalpy 

The AGN jet powers were calculated from the cavity en¬ 
thalpy and an estimate of the cavity system’s age. The age 


can be calculated in several ways (McNamara & Nulsen 
2007), but is generally assumed to be the buoyancy time¬ 


scale for the cavity to rise from the nucleus to its projected 
location. The power could be significantly higher if the cavi¬ 
ties are launched with a substantial initial velocity. Enthalpy 
is a more direct measurement. We assume the cavity fill¬ 
ing material is relativistic, so the ratio of specific heats is 
7 = 4/3, making the enthalpy 4pK ( Rafferty et al.|[2006 |. 
But it is unclear in the case of multiple cavity systems what 
the representative value to use would be. In Figurej^we plot 
the enthalpy of the cavity system nearest to the nucleus for 
each cluster versus the iron radius. This is a measure of the 
energy output of the AGN in the past 100 Myr or so. The 
dashed line is again the least-squares linear regression best¬ 
fitting. 


Rpe = (57± 30) X eSv (kpc). 


(3) 


where the enthalpy is in units of 10®® erg. The standard 
deviation is approximately 1.56 dex, somewhat larger than 
the power plot. However, as expected the relationship with 
increasing iron radius remains. Hydrodynamical simulations 


performed by Morsony et al. (20101 show the power, not the 


total energy, determines the volume an AGN can influence. 
Our result cannot confirm this conclusion, but jet power 
does seem to be more closely correlated. 

Finally, we consider how radio synchrotron power and 
iron radius compare to the mechanical jet power. Using prin¬ 
cipal component analysis, we can compare how different pa¬ 
rameters (jet power, enthalpy, bolometric and 1.4 GHz radio 
luminosity) correlate with the iron radius. For jet power, the 
correlation with iron radius is 0.92. In comparison, the cor¬ 
relation between enthalpy and iron radius is lower at 0.85. 
Finally, the lowest correlation with iron radius is the radio 
luminosity (bolometric and 1.4 GH), with a value of 0.74. 
Figurej^shows how much tighter the relationship is between 
iron radius and jet power than the relationship between ra¬ 
dio synchrotron power and jet power. The green crosses rep¬ 
resent the jet powers and the red diamonds and blue squares 
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Figure 4. - continued-. Top row: Hydra A and HCG 62. Middle row: A1795 and A1835. Bottom row: A2029 and A2052. 
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Figure 4. - continued: Top row: A2199 and Sersic 159. Bottom row: A2597 and M87. 


represent the bolometric radio luminosity and 1.4 GHz lu¬ 
minosity, respectively. This is to be expected to some degree 
because of the large scatter in the jet power vs radio syn¬ 
chrotron power correlation (Birzan et al.|[2004 Rafferty et 


al. 2006 Dunn & Fabian 20081. However, Figure shows 


that the jet power provides the best proxy for the average 
iron radius. 


4.4 Lifting Hot Gas by Radio Jets 

Large quantities of plasma are being lifted out of the centres 
of galaxy clusters, so a significant fraction of the bubble 
enthalpy must be expended lifting the gas. This quantity was 
originally calculated for Hydra A, and the work done lifting 
the plasma was found to be comparable to the enthalpy of 
the inner cavity pair (Kirkpatrick et al. |2009a (. We have 
carried out a similar analysis on the High Quality sample. 

We assume the outflowing gas is confined to a cylindri¬ 


cal volume with the same approximate radius as the cavity 
system. It is unrealistic to believe that the metals are uni¬ 
formly filling the conical regions shown in Figure used 
to define the average metallicity along the jets. Because we 
do not know the hlling factor, our best estimate is to as¬ 
sume at most the metals fill a volume with a width no 
greater than the observed cavity. We base this off of the 
highest resolution metallicity maps, such as M87 and Hydra 
A, where the plumes of metal-rich gas appear to be narrow. 
The length of the cylindrical volume is well known, extend¬ 
ing from the core to the iron radius. The width suffers the 
same uncertainties when measuring cavity sizes. This alone 
could further add a factor of 2 to 4 to the mass uncertain¬ 
ties. We further assume the off-jet profile is the undisturbed, 
azimuthally averaged prohle prior to the AGN outburst. Fi¬ 
nally, we assume the uplifted metals were originally pro¬ 
duced within the half-light radius of the BGG. 

The uplifted iron mass is calculated as the sum of the 
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Figure 5. Updated jet power vs. iron radius. This includes six 
more data points than the previous relationship presented in |Kirk-| 
Ipatrick et al.| | |2011[ |. The dashed line is the new best-fitting to 
the data. 


Figure 6. The enthalpy of the youngest cavity system for a clus¬ 
ter is plotted against the iron radius. The best-fitting to the data 
has slightly more scatter than that in Figure]^ 


iron mass within each radial bin of the outflow region, mi¬ 
nus the metallicity profile of the undisturbed azimuthally 
averaged profile. The uplifted iron mass is given as 


energy varies from a few percent to as much as 50%. Note, 
the iron mass measured here for Hydra A is smaller than we 
found earlier in Kirkpatrick et al. (2009a I owing to a more 
conservative estimate of the volume filled by the metals. 


Mpe = /Pe,© piVi{Zon,i — (4) 

i=l 


where p is the density of the gas assuming rie = 1.2nH, V 
is the volume of the outflow region, Zon and ZoS are the 
metallicity of the gas along the jet and orthogonal to the 
jet, and /pe,© is the iron mass fraction of the Sun. Assuming 
this iron came from the central region of the cluster, we can 
convert the iron mass values back into total uplifted mass. 
This is given as. 


Mgas 


Mpe 

Zofpe.Q ’ 


(5) 


where Zo is the central metallicity value. The energy that is 
required to uplift this gas can be estimated by calculating 
the difference in gravitational potential energy at the cluster 
centre and its current posiition. [Reynolds, Casper, fc Heinz| 
120081 gives the energy as 


AE = In 


( 6 ) 


where Cs is the sound speed, 7 = 5/3 is the ratio of specific 
heats for the X-ray emitting plasma, and pi and pf are the 
initial and final gas densities. The outflow masses, rate, and 
uplift energy for each cluster in the high quality sample can 
be found in Table These are very rough approximations 
due to the unknown filling factor of the estimated outflow¬ 
ing volume. The errors quoted here reflect the propagation 
of uncertainties due to the measured errors on density and 
metallicity. The total gas mass is roughly 3 orders of mag¬ 
nitude greater than the iron mass primarily due to the solar 
abundance fraction. Uplift energy as a fraction of the total 


4.5 Creating Broad Abundance Peaks Through 
AGN Activity 

AGN activity has been shown to be capable of redistributing 


gas in galaxy clusters (Simionescu et al. 2008 2009 Kirk¬ 


patrick et al. 2009a 2011 Gitti et al. 20111. Because this 


appears to be an efficient mechanism, there is reason to be¬ 
lieve AGN could in principle be responsible for the broad¬ 
ened abundance peak observed in cool-core clusters. We have 
tested the feasibility of AGN broadened abundance peaks by 
calculating the amount of energy required to move the entire 
mass of the abundance peak from the centre of the cluster 
to its current position using the same methods as described 
in the previous section. 

Azimuthally averaged abundance profiles have been cre¬ 
ated for the entire sample with 6 to 8 bins per profile, span¬ 
ning the field of view of the CCD. The large bin size allows 
for lower uncertainties and the spatial resolution to distin¬ 
guish the abundance peak from the background. The pro¬ 
files were fit using the method described in section [34] Only 
three clusters (2A 0335, A262, and A2052) were fitted us¬ 
ing a two-temperature model in the central region. The left 
panel of Figure is an ordinary example of the resulting 
profiles. Targets such as Perseus and M87, where the field 
of view only covers a small part of the cluster centre, have 
been excluded from this analysis. The dotted vertical lines 
indicate the effective radius of the BCG from the HyperLeda 
catalogu^ If material were not diffused outward, then the 


^ http://leda.univ-lyonl.fr 


© 0000 RAS, MNRAS 000, 000-000 


























































12 C. C. Kirkpatrick & B. R. McNamara 


Table 2. Jet properties 



MFe“ 

Mgaa® 

-^outflow 

F'uplift'* 

System 

(10^ Mq) 

(lOlo M©) 

(M© yr-i) 

( 10 ®® erg) 

A133 

9 74+0-75 
^•'^-0.69 

1-OU_o.25 

49± 21 

14-513 9 

A262 

f) 0951 +0.0125 

U.UOOi_Q 

0 01 94+0.0487 

U.U194_o 0421 

2.0 ±0.4 

0.0852tO;®279 

Perseus 

n 1 qq+0.016 

(J.ioy_Q 013 

0 154+0013 
U.104_o 010 

23 ± 12 

0.7801®;“®? 

2A 0335 


0 249+0-138 
0.^49_o.i45 

36 ± 14 

0.912l®-j“ 

A478 

2.98tl°l 

^•^O_o.78 

150 ± 62 

42 7+13.0 
^^+-14.5 

MS 0735 

9.40li;i° 

6 09 

0.uy_3 08 

150 ± 81 

243l?®| 

Hydra A 

0.893^0;!?® 

0 777+0.171 

0./ 1 <_0.157 

49 ± 19 

10.411“ 

HCG 62 

U.UZOD_0.0061 

n ni 99 + 0-0060 

u.uizz_o 0031 

3.8 ±0.4 

U.U^UZ_0.oo53 

A1795 

1.05l°J® 

0.793lgi™ 

77± 36 

6.64ll?“ 

A1835 


0 >^0 + 1-55 

0 . /9_i 40 

170 ± 86 

16ll®« 

A2029 

0.89lt0;«® 

0 472+0-220 
U.4^^_0 208 

49 ± 15 

IS 1+^’^ 
-‘-'^•-‘--4.6 

A2052 

n 1 91 +0.023 

O'.l2ii_o.022 

U.U040_o 0100 

2.2 ± 0.8 

0 974+0.084 

U.9/4_0.o82 

A2199 

0 061 0+0-0509 
U.UDlU_o.o501 

0 0*11 9+0-0261 
u.uoio_o 0257 

4.6 ± 1.9 

0 o-jQ+0.219 
U.9iy_o 211 

Sersic 159 

0 R94+0-277 

0l.O94_0.243 

0 501 +0-220 

U.OUl_o 193 

73 ±30 

9 Q7+1-04 

'^•9^_o.92 

A2597 

1 44+0-384 
J-.44_o 408 

1 04^+0.360 
-‘-•'^^-0.383 

98± 50 

14.511? 

M87 

n noc9+0.0039 
U.UOOZ_o 0039 

U.U2oy_o 0043 

0.92 ±0.18 

0-284lo®o®38 


“ Mass of uplifted iron 
^ Total mass of uplifted gas 
Outflow rate of the gas 

Energy required to raise gas from centre to its observed location 



Figure 7. Bolometric radio power(red), 1.4 GHz radio power 
(blue) and AGN jet power (green) have been plotted on the 
same scale. The corresponding coloured lines are the correlations 
best fit. The jet power, with a correlation coefficient of 0.92, is a 
stronger indicator of iron radius. The correlation coefficient with 
either radio power is 0.74. 


abundance peak should lie within this radius (De Grandi & 


Molendi||2001| |De Grandi et al.||2004| [Rebusco et al.||2005[ 

2006[[David fc Nulsen||2008 1. 


We first must dehne peak and background metallicity. 
We started by fitting a spline to each abundance profile. 


Taking the derivative of the spline ht, we searched for an 
inflection point in slope where the profile is transitioning to 
the background value. We choose this point as the turnover 
in slope closest to zero, where the slope is rapidly changing 
interior to this point and remaining near zero exterior to 
this point. We refer to this turnover point as the significant 
abundance peak radius (i?z)- The right column of Figure]^ 
is the change in slope plotted against radius. The dashed 
line represents Rz- 

To test whether our defined metallicity peak is more 
extended than the light, we have calculated the half-metal 
radius. The half-metal radius is the radius at which half the 
iron within Rz, by mass, is located. In Figurej^we present a 
histogram of the half-metal to half-light radius fraction. For 
the nine clusters where we have an effective radius measure¬ 
ment, we find that the half-metal radius is 1.5 to 3.5 time 
more extended than the light. If there was no broadening of 
the abundance peak, these values would be expected to be 
similar. 

With a peak radius defined, we can now estimate the 
approximate metallicity background. We use the average 
metallicity value of the profile points outside of Rz to be 
the subtracted background. The modified version of equa¬ 
tion [4] is now 

n 

Mfs = /pe,© PiVi{Zi — Zba.ck), (7) 

i=l 

where V is the volume of the spherical shells, Z is the total 
metallicity, and .Zback is the background metallicity. For the 
purpose of this analysis, the iron mass given here is assumed 
to have all originated from the BGG. Following the same as¬ 
sumptions we made for the mass uplifted by a single jet, the 
iron mass has been converted to a total gas mass using equa¬ 
tion]^ The final step is to calculate the energy to displace 
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Figure 8 . Left. The abundance profile of Abell 133. Profiles such as this were created for all clusters in the sample. The dotted line is 
the half-light radius of the BCG. The dashed line is the abundance peak radius Right. The derivative of the spline fit to the abundance 
profile. The abundance peak radius was derived here at the point where the slope approaches zero and remains mostly flat. 
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Figure 9. The fraction of half-metal radius to half-light radius. 
Half of the iron produced within the BCG is extended many times 
beyond its effective radius. The small subset of clusters used here 
are from the objects in Table 1 and Table 2 where effective radii 
measurements were available in the HyperLeda catalogue. 


the total mass of gas from the centre to its current observed 
location using equationj^ All values calculated here are pre¬ 
sented in Table The average background metallicity over 
the sample is 0.422 Zq and most clusters in the sample have 
displaced on the order of 10® Mq of iron. The energy to dis¬ 
place the iron ranges as low as 10®® erg for the galaxy group 
HCG 62 and as high as 10®® erg for massive clusters such as 
A1835 and MS 0735. 


In order to confirm that a single BCG can plausibly be 
responsible for the magnitude and width of a cluster abun¬ 
dance peak, we have determined the time-scales over which 
the material is created by SNe la, and how much the SMBHs 
must grow. If enrichment times are too long, this would in¬ 
dicate that some of the material is coming from another 
source (i.e., stripping), meaning a different mechanism be¬ 
sides AGN must be invoked to explain this phenomenon. 
Conversely, time-scales that are extremely short would sug¬ 
gest our abundance peak size calculations are grossly un¬ 
derestimated. To calculate enrichment time, the relationship 


used by Bohringer et al. (20041 takes into account two com¬ 


ponents, SNe la and stellar mass loss, that enrich the sur¬ 
rounding ICM with metals. The enrichment time to create 
all the metals in the abundance peak is given as 

1 .-1 Mpel/B0 

7Fe) (8) 


tenr = (lO'^^S’rjpe + 2.5 X 10“ 


where S = 0.15 SNU is the supernova rate (Cappellaro, 


Evans, & Turatto 19991, rjpe = 0.7 Mq is the SN la iron 
yield, ype = 2.8 x 10“'^ is stellar mass loss fraction of iron, 
and Lb is the total B-band luminosity referenced from Ta¬ 
ble Most clusters in the sample have little to no star for¬ 
mation in the BCG. An insignificant fraction of the total 
iron mass of the abundance peak is expected to originate 
from SN II ( Kirkpatrick et ^|2009a l. For most clusters in 
our sample, SN II only contributes at most 1% of the peak 
iron mass. However, A1835 stands out as exhibiting a high 
star formation rate between 100 and 180 Mq yr“^ (McNa- 
mara et al.||2006 1. At a star formation rate of over 100 Mq 
yr“^, the contribution of SN II would increase to over 10%. 
To account for this, equation]^ is modified to 


teni = Mpe[(10 ^^Sla^PeT 


Lb 


2.5x10 7Fe)7-h SFR7Pe,iiS'ii] , (9) 

Tbq 
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where SFR is the star formation rate, 7Fe,ii = 0.07 Mq is the 
SN II iron yield, and Sii = 0.01 SN Mg^. The enrichment 
time-scale for each cluster in our sample can be found in 
Table The enrichment times range from a few hundred 
Myr to a few Gyr. This is shorter than cluster ages so it is 
plausible for all the material to have been produced within 
the BCGs. 

To determine if enough energy is available, we measure 
the change in SMBH mass required to output the total en¬ 
ergy that went into uplifting all the metal-rich gas to form 
the abundance peak. The black hole grows by 

AMbh = (1 — e)Macc, (10) 


where e is the efficiency depending on the spin of the black 
hole and Macc is the accretion mass. The black hole spin is 
important here because it determines the radius of the last 
stable orbit, which in turn determines the binding energy 
that can be released. This efficiency can range from approx¬ 
imately e = 0.06 for a non-rotating Schwarzschild black hole 
to as much e = 0.4 for a maximal spinning Kerr black hole 
( King et al.|2002 I. For this analysis we have adopted a stan¬ 
dard value of e = 0.1. Accretion mass can be written in 
terms of the energy released, giving 


AMbh = (1 — e)—^ 


( 11 ) 


where E is the uplift energy, given from equation multi¬ 
plied by an additional efficiency factor. The extra efficiency 
factor is to account for the additional energy put into cavity 
expansion and shocks during an AGN induced outflow. We 
know this for each cluster by taking the uplift energy of a 
single outburst from Table and comparing it to the total 
energy output during the lifetime of the outflow. The aver¬ 
age age of outflows in our sample is approximately 10® yr 
and an average uplift-to-total energy fraction of 14%. The 
histogram of black hole mass growth are presented in Fig- 
ure[^and values reported in Table[^ Again, the errors here 
reflect the propagation of uncertainties due to the measured 
errors on density and metallicity. It is plausible that all of 
these SMBHs could have grown by this amount. The range 
of total black hole growth in our sample are consistent with 
previous estimations of SMBH masses ( Rafferty et al.|2006 1. 
The largest outburst in our sample, MS 0735, is on the same 


order growth of what McNamara et al. (20091 found to be 
the maximum growth the SMBH must undergo to create 
the observed cavity system. When considering the spin of a 
black hole, the values reported in Table[^would be reduced. 
Less accretion mass is needed in a spinning black hole to 
maintain the same output, resulting in less overall growth. 
At the maximal spin rate, the change in black hole mass for 
every cluster in our sample would be reduced by an order of 
magnitude. 

It may be surprising at first glance that MS0735, the 
system with the most energetic outburst, does not have the 
largest estimated SMBH growth. However, the mass and 
overall distribution of iron in a cluster core encodes the 
power output history of the SMBH. The integrated power of 
the AGN over time will generally exceed the instantaneous 
power. The estimated SMBH growth found for A1835 is the 
largest because its abundance peak has the largest iron mass 
and therefore the largest expended enthalpy. 

Radio AGN outflows may be responsible for the broadly 



10'’ 10® 10® 10^ 10® 10® io’° 

AMbh [ Msoi ] 


Figure 10. Histogram of the change in mass of the SMBH would 
have to take place to account for all the energy used in uplifting 
all the material in the abundance peak. 


distributed abundance peaks in cool core clusters. Across 
our sample, it takes on average 2 Gyr to produce all the 
iron enriching the abundance peaks through the process of 
SN la explosions, stellar mass loss, and core-collapse super¬ 
novae associated with star formation. The change in black 
hole mass associated with the energy required to uplift gas 
ranges from 10® - 10® M©. This range is less than the ex¬ 
pected masses of SMBH in large BCGs and is consistent with 
previous predictions of expected growth. Adding black hole 
spin to the calculation lowers the threshold for the amount 
of accretion material need for equivalent energy output. We 
believe because of these easily achievable standards, AGN 
are the prime mechanism for redistributing all the metal 
enriched gas associated with the BCG. 


4.6 Regulating Star Formation by Hot and Cold 
Molecnlar Outflows 


The hot outflow rates estimated here are comparable to or 


exceed the star formation rates found in BCGs (Rafferty et 
aL]|2006 1. Therefore, outflows of hot gas may play a role in 


offsetting cooling and regulating star formation. 

This conjecture is supported by recent ALMA discover¬ 
ies of what appears to be outflowing molecular gas in Abell 
1835 (included in this sample), and perhaps in Abell 1664 
(Russell et al.|2014 McNamara et al.|2014 1. Similar molecu¬ 


lar flows are present in the Perseus cluster (Lim, Ao, & Dinh- 
|V-Trung|2008[[Salome et al.|20lT| . The molecular gas flows 
are surprisingly large, with masses of 10® M© to lO'^® Mq. 
The hot and cold outflows may be physically related to each 
other. As the rising bubbles lift the metal-rich hot gas, the 
cooler ^ 1 keV gas may condense into molecular clouds as 
it rises outward behind the bubble (McNamara et ar]|2014 


|Li &: Bryan|2014| ). An alternative but more challenging pos¬ 
sibility would be that molecular clouds are dragged upward 
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Table 3. Properties of metal-enriched gas in galaxy clusters abundance peaks. 



Lb“ 

7 b 

^back 

E= 

Mpe'* 

AMbh'’ 

tenr"^ 

System 

10“ Lq 

Zq 

10®® erg 

lOl Mq 

10'^ Mq 

Gyr 

A133 

1.77 

0.467 

2331^^ 

52.4111 

Q qq+0.86 

0.cS0_2_78 

9 71 +0.12 
' -^-0.40 

A262 

0.417 

0.481 

3.95«;®? 

1.6311“ 

0 i42+^-®i^ 

U.035 

n qc;q+0-018 

U.900_o 039 

2A 0335 

2.56 

0.534 


7-86ll,“ 

1.6611“ 

0.28ll?:l?l 

A478 

4.06 

0.345 

185011260 

1081“ 

66.61“:" 

2.4311-1? 

MS 0735 

1.16 

0.346 


96.8l"l 

43.6l?l? 

7 63+^-!’^ 

+ 0'J_2.64 

Hydra A 

2.58 

0.345 



5.751111 

U.b5U_o 260 

HCG 62 

0.151 

0.333 

0.189l°;®“ 

0-188ll°o55 

u.UiU_o 007 

^•-LJ-^_0.033 

A1795 

1.69 

0.349 

11 70+^^® 

95.7l";°5 

42.011?? 

cr to + 1.06 

A1835 

3.16 

0.324 

57381“®° 

1721^1 

20611?? 

4.9811:?? 

A2029 

1.64 

0.455 

5741^6 

66.2l“+ 

2o.6ll?:l 

3.6911:11 

A2052 

0.773 

0.526 

69.5l«;® 

17.8111 

^•OU_o 94 

91-1 +0.54 
^•-*--*--0.41 

A2199 

2.49 

0.492 

47.ll“i® 

11 7+^-^ 

^ -0.2 

1 69"^^'^^ 

-‘-•o^_0.08 

0.4301?-°?? 

Sersic 159 

1.34 

0.456 

36.4l|° 

10 1+^-^ 

1 “^1 +^' 1^ 
J--oi_o 10 

0.6861?:??? 

A2597 

1.57 

0.439 

1031° 

15 4+0.4 

q "70 + 0.21 
‘ ^-0.20 

0.8961?:??? 

A85 

1.98 

0.470 

221+_lt 

41 1+^-^ 
^-*-•-*--1.8 

7 94+^-®^ 

< •y^_0.78 

-L"»'^_0.08 

PKS 0745 

0.254 

0.494 

22A+11 

o 9 c:-I-0.47 
8.25_o,45 

0.805l?:ll? 

2 97+0-17 
-0.16 

4C 55.16 

1.12 
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386lg7 

71.411:1 

i3.9i?:i 

c- oq-l-0.49 
5.83_(, 45 

RBS 797 

— 

0.351 

1660+"° 

1221" 

59.6111? 
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A2390 

— 

0.366 

2511“ 

61.8111 
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1.24 

0.418 

ou.i_i8 1 

56.3l« 

2.1611?? 

4 i»+0-34 
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2.07 
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3 10+0.27 
'^•-‘-'^-0.25 
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— 

0.416 


80.4l?'l 
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— 

MKW 3S 

1.07 
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1501" 

31.811:1 
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o.oy_o 51 

9 yq + U.lC5 

‘'^-0.13 
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56211“ 

69.9111 

28.31?:? 

10 2"^^’^ 

Cygnus A 
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61 

- 17.7 

9 93+5-87 

y.yo_5 75 

6 71 + 1-11 
' 1--1.92 
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2.21 

0.607 

40.411° 

14 7+^-’^ 
l^+-0.6 

^■^^-0.23 

u.Duy_o 027 


“ S-band luminosity derived from HyperLeda catalogue 
^ Background metallicity 

Estimated uplift energy -|- cavity/shock energy 
Abundance peak iron mass 

® Change in black hole mass needed to account for total output of energy 
■f Enrichment time to create the observed iron mass 


in the wakes of the more massive hot outflows (McNamara 
et al.|20l4 1 . The outflowing hot and cold gas, perhaps asso¬ 
ciated with a circulating flow ([Mathews et al.||2003| |Math-| 


ews, Brighenti, fc Buote]|2004|]McCarthy et al. 20081 may 

be delaying or removing gas otherwise destined to fuel star 
formation in the BCG. Figure [TT] shows that the total out¬ 
flowing gas masses lie between ICr Mq to 10^^ Mq across our 
sample. We have calculated the associated flow rates for each 
cluster assuming the gas is traveling at the bubble buoyancy 
speeds but no greater than the sound speed of the ICM. The 
hot outflow rates lie between a few solar masses per year for 
jet powers of ~ 10^® erg s“^ to more than 100 Mq yr“^ for 
jet powers exceeding ~ 10^® erg s“^. A crude relationship 
between jet power and outflow rate is shown in Figure [T2| 
The dashed line indicates the least-squares linear regression 
of the form, 

Moutflow = (22 ± 16) X (Mq yr-^). (12) 


Here, jet power is in units of lO"^'^ erg s“^. The scatter about 
the mean is large, with a standard deviation of approxi¬ 
mately 1.47 dex. Comparing the outflow rate to the enthalpy 
doesn’t not improve the scatter. 


We compare the outflow rates to the nominal cooling 
rates of the gas calculated as. 


• _ 2 /rrWpLx 

“°‘ 5 fcbTx ’ 


(13) 


where Lx is the cooling luminosity from [Rafferty et al.j 
( [2006[ ) and Lx is the average gas temperature within the 
cooling radius. In Figure we present a histogram of the 
ratios of outflow rates to the cooling rate. While the un¬ 
certainties are admittedly large (as high as 50%), the figure 
shows that outflows are able to remove only 10% to 20% on 
average of the gas that is expected to cool. The higher frac¬ 
tions in the histogram are not associated with with higher 
jet powers, and likely reflect uncertainty in the outflowing 
gas masses. This figure shows that hot outflows alone would 
be unable to remove all of the gas that is expected to cool 
and to accrete onto the central galaxy. Instead, they may 
be one part of many processes which regulate the rate of 
cooling gas which fuel star formation and the nuclear black 
hole. 
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l^outflow [ I^sol ] 


Figure 11. Distribution in total mass outflowing from each clus¬ 
ter in the high quality sample. 



^outflow [l^sol ] 

Figure 12. Jet power vs. mass flow rate for systems in the high 
quality sample. The dashed line is the best-fitting to the data. 


5 SUMMARY 

The 16 clusters of the High Quality sample we have pre¬ 
sented here reveal higher metallicities along projected tra¬ 
jectories of their X-ray cavities. We interpret this effect as 
the signature of massive, hot outflows of metal-enriched gas 
produced within the BCG. At the same time the metallic- 
ity maps reveal a great deal of complexity, indicating that 
other processes, including returning flows of gas, may be 
taking place. The outer iron radius is found to scale with jet 
power as Pjir S'lid with cavity enthalpy as The jet 

power relation shows slightly lower scatter. 



I^outflow ! I^cool 


Figure 13. The fraction of outflowing material to the expected 
total cooling rate of the cluster. The redistribution of material by 
AGN alone is insufficient to balance a cooling flow. 


The mean hot outflow rates are typically tens of solar 
masses per year and upward of 100 M© yr“^ in the extreme. 
These flow rates are generally insufficient to balance the 
radiative cooling rate of a cluster. Recirculation of cooling 
gas accounts for only 10% to 20% of the radiative cooling 
rate, so that heating must be more important. Nevertheless, 
the outflow rates are comparable to star formation rates in 
BCGs, and may have a significant impact on fueling star 
formation and the AGN itself. 

Using the entire sample, we found that over the life 
cycle of AGN activity, hot outflows may alone be responsible 
for the broadening of abundance peaks in cool-core clusters, 
and effectively transport metal-enriched gas from the BCG 
to great distances in cluster atmospheres. On average, only 
~ 15% of the energy released by the AGN is required to lift 
the gas from the centre to the observed altitudes. Standard 
conversions between free energy and rest mass indicate that 
the SMBHs masses must have increased between 10^ and 
10® Mq while driving the outflows. This mass increase lies 
comfortably within the estimated SMBH masses of BCGs, 
and is therefore a plausible explanation for the extended 
abundance peaks in clusters. 
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